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\Q ; Various types of galaxies observed in the cosmological scales show P-Cygni type pro- 

files in the Lyct emission lines. The main underlying mechanism for the profile formation 
in these systems is thought to be the frequency redistribution of the line photons in an 
expanding H I component surrounding the emission source. 

A Monte Carlo code is developed to investigate the Lyct line transfer in an optically 
thick and moving medium with a careful consideration of the scattering in the damping 
wing. Typical column densities and the expansion velocities investigated in this study are 
N HI ~ 10 17 " 20 cm" 2 and AV ~ 100 km s _1 . 

The main features in emergent line profiles include a primary emission peak in the red 
Q-f part and a much weaker secondary peak in the blue part. We investigate the dependence of 
the profile on the kinematics and on the column density. The primary peak recedes to the 
red and the width of the feature increases as Nhi increases. It is also found that symmetric 
double peak profiles are obtained in the static limit, where bulk motion is negligible. 

The P-Cygni type profile in the Lya emission line of DLA 2233+131 (z = 3.15) is 
noted to be similar to those found in the nearby galaxies and distant (3 < z < 3.5) H II 
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galaxies. Our numerical results are applied to show that the DLA system may possess an 
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expanding H I supershell with bulk flow of ~ 200 km s _1 and that the H I column density 
Nhi is approximately 10 20 cm -2 . 

From the observed Lya flux and adopting a typical size of the emission region ~ 1 kpc, 
we estimate the electron density of the H II region to be ~ 1 cm -3 and the mass of H II 
region ~ 10 8 Mq. We also conclude that it requires ~ 10 3 05 stars for photoionization, 
which is comparable to first-ranked H II regions found in nearby spiral and irregular 
galaxies. 

We briefly review the physical quantities of several astronomical objects which may 
possess similar kinematic properties to those of DLA 2233+131. We point out the fact that 
this kind of outflowing media are often accompanied by the sites associated with active 
star formation. 

Key words : Radiative transfer - :Monte Carlo - QSO:damped Lyman alpha sys- 
tem:individual:DLA 2233+131 - cosmology:galaxy formation 
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1. Introduction 



A typical high-redshift quasar spectrum shows a variety of absorption features blue- 
ward of the broad Lyct emission line, most of which are not associated with the quasar itself. 
These absorption systems are thought to be related with the neutral hydrogen and metal 
elements in the intervening medium along the line of sight to the quasar, and therefore 
constitute an ideal subject for studying the cosmological structure formation processes. 
Depending on the neutral hydrogen column density Nhi, the absorption systems are usu- 
ally classified into three classes, i.e. damped Lyct (hereafter DLA) systems, Lyman limit 
systems, and Lyct forest systems. 

With the advent of the Hubble Space Telescope and 10 m class telescopes such as 
Keck I, II, detailed analyses of absorption line profiles have been performed, and provided 
very important kinematic information about the absorption systems. 

Recently Prochaska & Wolfe (1997) analyzed the absorption line profiles of the 17 
DLAs obtained with the HIRES eschelle spectrograph on the Keck telescope and concluded 
that the kinematics of the DLAs supports the hypothesis that they are rapidly rotating 
"cold" disks (see also Haehnelt et al. 1997). 

There have been a number of DLA galaxy candidates possibly exhibiting Lya emission 
features in their spectra (e.g. Hunstead et al. 1990, Moller & Warren 1993, Pettini et 
al. 1995 etc. ). However, the first clear Lya emission feature was detected in the spectra 
of the DLA 2233+131 (Djorgovski et al. 1996). 

Djorgovski et al. (1996) reported an isolated galaxy (or a protogalaxy) 2.3" (or 17 kpc) 
away from the line of sight to the quasar, QSO 2233+131 by detecting the Lya emission 
having the same redshift as DLA 2233+131. They also noted that the emission profile is 
asymmetric to the red and that there is a possible absorption feature to the blue. The 
width of the absorption feature is measured to be ~ 200 km s _1 . It is not certain that the 
absorption is related to a very near Lya forest, as was pointed out by several researchers 
(Djorgovski et al. 1996, Lu et al. 1997). However, the similarity of the profile to the P- 
Cygni profile leads to an interesting possibility that an expanding medium envelops the 
emission source, which can be thought to be a supershell seen in nearby galaxies (see section 
4. of Irwin (1995) for a review of outflowing structures). Since the first firm identification 
of DLA candidate DLA 2233+131, Djorgovski et al. (1997) found several new candidates. 

In addition to this direction, there are at least two other groups adopting efficient 
programs to find out primeval galaxies. 

Firstly, Giavalisco et al. (1996) found that there are a number of star- forming galaxies 
at 3.0 < z < 3.5 possessing compact spheroidal cores often surrounded by low surface 
brightness nebulosities. Morphologically similar objects are also found in the Hubble Deep 
Field (Steidel et al. 1996a). Moreover, Steidel et al. (1996b) argue that the spectra of 
galaxies with z > 3 are remarkably similar to those of nearby star-forming galaxies and 
often show Lya lines with P-Cygni profiles characteristic of expanding H I media. The H I 
column density is inferred to be 10 20 cm -2 in these systems. 

One excellent example of nearby star-forming galaxies is provided by Legrand et 
al. (1997) and Lequeux (1995) who performed a detailed analysis of the emission pro- 
files of the star- forming galaxy Haro 2. They suggest that these P-Cygni type profiles are 
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formed by the expanding H I shell surrounding an H II region. At least three more galaxies 
showing such profiles are found in the literature (Kunth et al. 1996). 

Recently Steidel et al. (1997) used the Lyman break criterion to obtain ~ 80 new 
candidates that are similar to those found previously. What is interesting is that ~ 75% 
of them exhibit P-Cygni type Lyct profiles (Pettini et al. 1997). Lowenthal et al. (1997) 
also observed the spectra of 24 remote galaxies selected from the Hubble Deep Field. 

Secondly, using the flux magnification via gravitational lensing by foreground clusters 
of galaxies, very faint and remote galaxies can be found. Trager et al. (1997) found a 
multiply-imaged galaxy at z ~ 4 gravitationally lensed by the rich cluster CL0939+4713, 
which are the same kind of those discovered by Lowenthal et al. (1997). Using the same 
tactics, Franx et al. (1997) discovered the currently remotest galaxy at z = 4.92. Recently 
Frye et al. (1997) found the remote galaxy at z = 4.04 which is lensed by the cluster of 
galaxies, A2390 (z = 0.23). Surprisingly, ~ 75% of the primeval galaxies described above 
exhibit P-Cygni type Lya emission lines, which implies the ubiquity of an expanding HI 
envelope around HII region. 

So far, many investigations about the P-Cygni profiles have been concentrated on the 
outflowing systems possessing rather moderate optical depth. The Sobolev approxima- 
tion is regarded as a powerful method to understand the behavior of resonantly scattered 
photons in an expanding medium with a bulk flow much larger than the thermal speed 
(Sobolev 1960, Mihalas 1978). However, when the medium has high optical thickness, the 
scattering in the damping wing is not negligible any more and the validity of the Sobolev 
approximation becomes questionable. 

In this paper, we develop a Monte Carlo code for the radiative transfer of Lya photons 
in an optically thick and expanding medium that is expected in the H II superbubbles in 
star-forming galaxies. Extrapolating the kinematical information from our line profile 
study, we try to find a possible relation of the P-Cygni type absorption profile to the 
existence of superbubble in the primeval starburst galaxies. 

The paper is composed as follows. In section 2 we describe the Monte Carlo method 
to deal with the Lya transfer in a thick moving medium. In section 3 we present the 
main results. We apply our results to the absorption features in the Lyct emission of DLA 
2233+131 in the following section. Discussions and the implications for the identity of 
DLA 2233+131 are given in the final section. 

2. Monte Carlo Simulations 



2.1 Model Description 

In this subsections we give a brief model description and a basic atomic physics for 
the Monte Carlo code, which computes the profile of the emergent Lyct photons scattered 
inside an optically thick and expanding medium. 

The photon source is assumed to be located at the center of the spherical scattering 
medium which is truncated at r = r ma x- We assume that the constant density distribution 
n{r) = no throughout the scattering region. We only consider an isotropic and point-like 
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line emission source in this work, which is not a bad approximation for emission sources 
such as H II regions surrounded by H I regions. 

The velocity field of the scattering medium is assumed to be given by a Hubble-type 



flow, i.e. , 



v = H v r, (1) 



where H v is the bulk flow velocity gradient, v is the bulk velocity, and r is the distance 
from the center of the velocity field. In this type of velocity field, the distance is more 
conveniently measured by a parameter s defined by 

s = H v r/v th , (2) 

where v t h is the thermal velocity. 

If we introduce the Doppler width Aujy defined by 

Au D = u — , (3) 

c 

then the frequency of a line photon is also conveniently described by the normalized fre- 
quency shift x defined by 

x={y- vq)/Av d . (4) 

Because the scattering medium is assumed to be optically thick, scattering in the 
damping wing should be considered with much caution. Therefore we introduce the damp- 
ing coefficient, a, normalized by Ai/£>, i.e. , 

a = T/4tvAu d , (5) 

where V is the damping constant. The scattering cross section in the rest frame of the 
scatterer is then given by 

yre 2 T/47T 2 
a u = a{x) = Josc- 



m e c JObC (v-v ) 2 + (r/4n)2 
ire 2 la 
m e c tyAud x z + cr 

where f osc is the oscillator strength, m e is the electron mass (Rybicki & Lightman 1979). 

The scattering optical depth Ti2(x) of a line photon with frequency shift x and wave 
vector ki corresponding to the distance from the position ri to r2 is computed to determine 
the subsequent scattering position in the Monte Carlo calculation. Due to the Hubble-type 
velocity field and a constant density field inside the scattering medium, we get the same 
Hubble-type flow if we make a translation so that ri coincides with the new origin. The 
line center frequency v'o{r) in the observer's frame differs from that in the scatterer's rest 
frame by the Doppler shift, which is given by 

v'o( r ) = ^o(l - H v r/c + vioc/c), (7) 
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where r is the distance from the new origin ri and v\ oc is the local velocity of the scatterer 
due to the thermal motion. 

The normalized frequency shift x'(r) is now simply given by 

x'(r) = x — s + u, (8) 

where u = vi oc /v t h- Therefore the scattering cross section a(x) becomes a function of the 
position due to the coupling of the line center frequency and the position determined by 
the expansion of the medium. 

The thermal distribution of the scatterers is assumed to be given by a Maxwellian 
distribution, that is, 

n = n Q J dvi oc f(yioc), (9) 

i f VLoc ) 2 i _ 2 
where f{vi oc ) = — e v v th = —=± — e~ u is the Maxwellian distribution function. 

Combining Eqs. 1-9, we obtain 



/CXD PS 
du \ ds' n f(u)a{x'(s')] 
-oo Jo 

= tq du e~ u 



, , U + X \ _■< / u + x — s 

tan ( I — tan 1 

a 



(10) 



Here, to is the Sobolev-type optical depth defined by 

yre 2 n a A , . 

(Sobolev 1960). 

The inverse transformation of Eq. 10 is obtained to find the normalized path length 
s corresponding to a given scattering optical depth t\i. The escape optical depth r esc is 
also obtained by setting s = s max , where 

Smax = H v r max /V t h = AVbulk/Vth- (12) 

Here, AVbuik is the difference of the bulk velocities at the center and the boundary of the 
scattering medium. 

The conversion of the H I column density and the Sobolev type optical depth r is 
then given by 

_ f Nhi \ / H v ■ r max \ / f osc \ 

T °~ U.l x 10 18 cm" 2 M 100 km s" 1 / V 0.4162 ' K ' 



2.2 Monte Carlo Approach 
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In this subsection, we present a detailed description of the Monte Carlo procedure. 
There are a few approaches to the resonance line transfer in an optically thick and static 
medium (Osterbrock 1962, Adams 1972, Harrington 1973, Sengupta 1994, Gould & Wein- 
berg 1996). However, there are also studies on the line formation in a thick expanding 
medium (e.g. Rybicki & Hummer 1978). 

The Sobolev approximation has been one of the most favored methods in dealing with 
the line transfer in an expanding medium. However, the validity of the Sobolev method is 
limited to the cases where the bulk flow is much larger than the thermal velocity and the 
optical depth of the medium does not greatly exceed unit. This is because the scattering 
in the damping wing is no more negligible in an optically thick medium where the typical 
line center optical depth r c > 10 4 . Therefore we develop a Monte Carlo code in which the 
scattering in the damping wing is carefully considered in a thick moving medium. 

The Monte Carlo code begins with the choice of the frequency and the propagation 
direction of the incident photon from an assumed Lyct profile. 

Then we determine the next scattering site separated from the initial point by the 
normalized propagation length s defined by Eq. 2, which corresponds to the optical depth 
of t = — ln(R), where R is a uniform random number in the interval (0, 1). Due to the 
unwieldiness of the inverse relation of Eq. 10, we tabulate the normalized path length 
s(x, r) in advance and look up the table to find s by the linear interpolation. For our 
numerical computation, the range of x is taken to be from — 4s max to 4s max with a step of 
Ax = 0.5, and r runs from to As max with a step of 0.002, where s max = 10 is assumed. 

In Fig. 1 we present the s(x, r) table as a surface plot. We restrict the propagating 
length s to be < s < 2s max , and set s = 2s max + 1 if the true value of s(x, r) exceeds 
2s max . A U-shaped distribution of s is seen in the figure, which is naturally explained by 
the fact that the scattering medium is transparent for extreme blue and red photons. 

The emitted photon traverses a distance s found by the above procedure and is scat- 
tered off if s < 2s max . In this scattering event the frequencies of the absorbed photon 
and the re-emitted one in the rest frame of the scatterer should be matched. For the case 
of an optically thin and expanding medium, this frequency matching is combined with 
the Sobolev approximation to yield the absorption profile in the form of the Dirac delta 
function. 

As mentioned earlier, in contrast with the case of a thin medium, for a very thick 
medium the scattering in the damping wing is not negligible. A good care need to be 
exercised to distinguish the scattering in the damping wing from the resonance scattering, 
because they show quite different behaviors in the properties including the scattering phase 
function and the polarization (Lee 1997, Lee & Blandford 1997, Blandford & Lee 1997). 

Because the natural line width is much smaller than the Doppler width, the local 
velocity of the scatterers that can resonantly scatter the incident photon is practically a 
single value. However, when the scattering occurs in the damping wing, the local velocity 
of the scatterer may run a rather large range. Therefore, in order to enhance the efficiency 
of the Monte Carlo method, it is desirable to determine the scattering type before we 
determine the local velocity u of the scatterer. 
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We present a more quantitative argument about the preceding remarks. Under the 
condition that a given photon is scattered off by an atom located at a position s, the local 
velocity component u along the direction kj is simply given by 



f( U ) K 7 M2 i 2 ' 

(it — o) 2 + a 2 

where 6 = s — x. The normalization condition is used to get 



(14) 



/(«) = 



7T 



-#(a, 6) 
a 



(u - 6) 2 + a 2 

Here, the Voigt function H(a, b) is evaluated by a series expansion in a, i.e. 

if (a, 6) = # (a, b) + aH^a, b) + a 2 H 2 (a, b) + a 3 H 3 (a, b) H , 



(15) 



(16) 



where H n (a, 6), n = 0, 1, 2, 3 are tabulated by Gray (1992). 

Because of the smallness of a, the function / has a sharp peak around u ~ b, for which 
the scattering is resonant. Therefore, the probability P r that a given scattering is resonant 
is approximately given by 



/oo 
d(Au) 
-oo 



(Aw) 2 + a 2 la 



n 



-H(a, b) 



H(a,b)' 



(17) 



The probability that a scattering occurs in the damping wing is 



P nr = 1 - P r . (18) 

In the code we determine the scattering type in accordance with the scattering type 
probabilities P r and P nr . If a scattering is chosen to be resonant, then we set u = b. 
Otherwise, the scattering occurs in the damping wing, and u is chosen according to the 
velocity probability distribution given by Eq. 15. 

We give the propagation direction k^ of a scattered photon in accordance with the 
isotropic phase function. The scattered velocity component v± perpendicular to the initial 
direction k^ on the plane spanned by kj and k/ is also governed by a Maxwell-Boltzman ve- 
locity distribution, which is numerically obtained using the subroutine gasdevQ suggested 
by Press et al. (1989). The contribution Ax of the perpendicular velocity component v± 
to the frequency shift along the direction of kf is obviously 

Ax = v ± [l-(k l -k f ) 2 ] 1 / 2 /c. (19) 



Therefore, the frequency shift Xf of the scattered photon is given by 

x f = Xi-u + u(k t ■ k f ) +v±[l- (ki • kf) 2 ] 1 / 2 , (20) 



where xi is the frequency shift of the incident photon. 
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In each scattering event the position of the scattered photon is checked and if it is out 
of the medium we collect this photon according to its frequency and escaping direction. 
The whole procedure is repeated to collect typically about 10 3 photons in each frequency 
bin. 

3. Results 

3.1 Escape Probability 

In Fig. 2 we show the emergent Lya profile from a thick expanding medium. The 
horizontal axis represents AA/AA^ = —x, and the vertical axis stands for the relative 
flux. Note that the wavelength devitaion is the negative of the frequency deviation and 
AXo is the Doppler wavelength. The emission source lying in the center of the scattering 
medium is assumed to be given by a Gaussian profile proportional to ^-^ x l' la ^) ; where the 
width o x is set to be 5. The use of a Gaussian profile may be justified by the observed 
Gaussian line profiles for other emission lines such as Ha or [O III] (Forbes et al. 1996, 
Legrand et al. 1997). We choose v t h = 10 km s -1 , s ma x = 10, which corresponds to 
AV6 u /fc = 100 km s -1 . We also choose r = 10 4 that leads to N HI = 4.1 x 10 18 cm -2 given 
the above choice of s max . 

The main features in the emergent profile in Fig. 2 include a primary peak in the 
red part and a much weaker secondary peak in the blue part. One may understand 
this behavior qualitatively by considering the escape probability P eS c, with which a given 
photon escapes from the region to an observer. We note that the single-scattering escape 
probability is given by 

Pesc = e- T ", (21) 

where the escaping optical depth r x is obtained from Eq. 10 setting s = s max . 

In Fig. 3 we show the single-scattering escaping optical depth and the escape prob- 
ability. The dashed line represents t x , and the P esc is shown by the dotted line. The 
dot-dashed line is the input emission profile at the origin given by a Gaussian and the 
solid line is the product of the initial emission at the origin and the escape probability. 
The analysis has a simple interpretation of the emergent profile without being scattered 
from the source to the observer. 

Because the medium is not optically thin, the escape probability analysis is a very 
poor approximation to the emergent profile which is severely affected by photons scattered 
a large number of times. However, the qualitative nature of the red asymmetric profile 
is obvious from the figure. The Sobolev optical depth tq ^> 1 for the incident photons 
with frequency shift x is satisfied in —3 < x < s max + 3. Therefore, we expect a severe 
discrepancy in this regime. 

In Fig. 4 we show the emergent profiles for the case where the incident profile from 
the emission source is given by a flat continuum. We use the same parameters as in Fig. 2. 
The broad absorption feature is formed in the frequency range of — 3 < x < s max + 3 with 
a sharp peak on the red side. 

3.2 Line Formation Mechanism 
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When the static medium has a moderate line center optical depth r c < 10 5 , the 
radiative transfer is dominated by the diffusion in frequency space, and the spatial diffusion 
plays a rather minor role. Hence each photon is scattered many times in the vicinity of 
its source and gets a sufficient frequency shift just before it escapes, which Adams (1972) 
described as "a single longest flight." For a thicker medium with r c > 10 5 , the damping 
wing scattering becomes important, and the photons emerge via the spatial diffusion. This 
process was called "a single longest excursion" by Adams (1972). He investigated these 
processes using a Monte Carlo method and showed that the mean scattering number < N > 
before escape is given by 

N ir c V^nF c forr c <10j (22) 
\ I^tt 1 / 2 ^ for r c > 10 5 . V ; 

In the case of an expanding medium, using the Sobolev theory in a resonance shell 
having a width of velocity difference of order vth, the effective optical depth is the Sobolev 
scattering optical depth tq that is the total scattering optical depth in the limit AV — > 0. 
However, this hypothesis works only when ro < 10 3 , in which case a frequency diffusion 
into the red wing leads to a direct escape from the scattering medium without further 
scattering in the damping wing. 

When To > 10 4 , the damping wing scattering optical depth is not negligible and the 
line photon can be captured with significant probability despite the frequency diffusion 
into the red wing. The subsequent scatterings tend to reduce the frequency deviation of 
the line photon as the restoring force in the frequency space operates (see Adams 1972). 
Therefore, the line photon is again resonantly scattered after several wing scatterings and 
severe spatial movement. In this case, the total scattering number is dominantly affected 
by the total column density of the scattering medium, which is proportional to the product 
of ro and s max according to Eq. 13. 

Therefore, when ro = 10 4 and s max = 10, the scattering number before escape is 
expected to be of order 10 5 by regarding the resonance shell as a static medium with a 
line center optical depth r c ~ ro (Adams 1972). It is checked that the scattering number 
increases almost linearly with TQS max . 

In Fig. 5a is shown a typical scattering behavior of a photon with x ~ — 5. The hori- 
zontal axis represents the scattering number. In the bottom panel is shown the frequency 
variation of the photon. In the middle panel, we show the radial distance of the photon 
from the center of the sphere with radius s = 10. The scattering type is shown in the top 
panel in which we denote a wing scattering by and a resonance scattering by 1. 

Near the 360th scattering, the frequency of the photon deviates much from the center 
frequency, and several scatterings in the damping wing follow. During this process, the 
photon propagates a large distance in real space and becomes resonant again by the effect 
of the restoring force (in frequency space) stated by Adams (1972). When the photon 
escapes from the medium, it experiences several wing scatterings, and the frequency of the 
photon falls on the red part. 

In Fig. 5b we show a typical scattering behavior of a blue photon with x ~ 5. Initially 
it experiences several damping wing scatterings and propagates easily with a large free 
path until it hits the resonance zone. Once it enters the resonance regime, the photon 
experiences the similar processes to those of the red photons described above. 
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On the other hand, majority of the extremely blue and the extremely red photons 
escape the medium suffering from at most several damping wing scatterings. 

3.3 Dependence on Kinematics and Column Density 

In this subsection we investigate the dependence of the emergent Lya profile on the 
kinematics and the column density of the scattering medium. Our numerical results are 
summarized in Figs. 6 and 7a,b. 

The kinematics of the scattering medium is characterized by the parameter s max that 
measures the bulk velocity scale. In Fig. 6 we show the emergent profiles for s max — 
0.1, 1, 2, 5, 10 with the Sobolev optical depth r = 10 4 and the width of the initial profile 
o~ x = 5. Here all the profiles have the same number of incident photons. It is particularly 
notable that the slope at the blue edge of the primary peak does not change as s max , 
which strongly implies that this quantity can be very useful to put a strong constraint on 
t . This point is discussed in more detail in the following section. 

As Smax gets smaller, the strength of the secondary peak approaches that of the 
primary peak. In the limit of s max — > 0, we obtain a symmetric double peak profile, 
and this is in agreement with the result obtained for the case of a static medium by 
Adams(1972). A very turbulent medium or a static medium may provide a negligible 

Smax- 

In Fig. 7a, we show the emergent profiles for moderate column densities, To = 10 2 , 10 3 
and 10 4 . Here, we fix s max = 10 and use an initial Gaussian profile having a width of 
o~ x = 5. Because all the profiles have the same number of incident photons, the profiles are 
normalized to have the same area. The thick solid line indicates the emergent Lya profile 
for to = 10 4 which corresponds to Nhi = 4.1 x 10 18 cm -2 , the dotted line for To = 10 3 , 
and the dashed line for r = 10 2 . 

The shapes of emergent profiles show a systematic change with the Sobolev optical 
depth. The primary red peak recedes to the red, as the column density increases. And the 
height of the primary peak decreases and the width increases. As To increases, the range of 
frequency diffusion increases and therefore the width of the primary red feature increases. 

In the cases of high optical depths r > 10 4 , the profile of the emergent flux changes 
in a similar way. In Fig. 7b we show the results for the cases tq = 10 4 , 10 5 , and 10 6 . Up 
to To = 10 5 the blue edges of the red primary peaks appear to occur at the same location 
x ~ — 2. The slope at the blue edge of the primary red peak decreases as the Sobolev 
optical depth, that is, the overall profile changes from a sharp primary feature to a broad 
hump as r increases. In contrast with the cases of Fig. 7a, the secondary blue peaks show 
remarkable changes. As To increases, the double peak profile becomes more symmetric, 
which is obtained in the static limit. This shows that the kinematics plays a minor role as 
the optical depth increases. 

For an extremely thick medium, To > 10 6 , the feature appears to have two broad humps 
outside the absorption trough. The strength of the red part is now almost similar to that 
of the blue part, which highly implies that the anisotropy introduced by the kinematics 
has little effect on the line formation. The enhanced width of the feature can be a good 
measure of the column density of the scattering medium. 

It is noted that the strength ratio of the primary red part to the secondary blue 
part is observationally important and that several sets of (to, s max ) may give the same 
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strength ratio. In the following section we give a further discussion about observational 
implications. 

3.4 Observational Implications 

We are primarily concerned with the quantities that can be observationally measured 
and put meaningful and strong constraints on tq and s max that specify the column density 
and the kinematics of the expanding H I medium. 

We propose that the slope at the blue edge of the primary part and the location of 
the primary peak can be one important indicator of To under the assumption that the 
spectroscopy can be performed with enough resolution and signal to noise ratio. 

In our work, the bin size for collecting the emergent photons is Ax = 1, that is, the 
bin size corresponds to the thermal velocity vth = 10T^ 2 km s _1 . The normalization is 
done by setting the total number of line photons to be 10000. Therefore, the slope at the 
blue edge is defined by the ratio of the photon number collected in the bin corresponding 
to the maximum peak point to the number of bins between the location of the blue edge 
and the maximum peak point of the primary red peak. 

As shown in Fig. 7b, for r = 10 4 , s max = 10, the slope at the blue edge of the 
primary part is given by 

3800 photons , . 

1 = TbTn ' [2S) 
and similarly for ro = 10 5 it is given by 

470 photons . 

' = 1 an ■ (24) 

It is particularly noted that these slopes are meaningful only when the total number of 
Lya photons throughout the feature is normalized to 10000. It is also interesting that 

I < 2 °P hotons . (25) 
1 bin 

for r = 10 6 . 

On the other hand, s max may be constrained by the quantities including the peak- 
to-peak distance, or equivalently the width of the absorption trough and the ratio of the 
strengths of the primary peak and the secondary peak. However, a good caution is needed 
because this quantity can be a poor indicator when a huge column density is involved. 

With lower quality spectroscopy we may use equivalent widths to characterize the 
emission profile in a similar manner which Ringwald & Naylor (1997) used in their analy- 
sis of the lines from the cataclysmic variable, BZ Camelopardalis. If we measure observa- 
tionally the equivalent widths EWi, EW2, and EW a b s corresponding to the primary red 
part, secondary blue part and the absorption trough respectively, then the profile fitting 
procedure will be reduced to minimize the differences of the equivalent widths obtained 
from numerical simulations. 



4. Application to DLA 2233+131 
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DLAs are believed to be associated with a disk component of a normal spiral galaxy 
or a dwarf galaxy from the considerations of the similarity of the H I column density 
to a typical galactic value and the width of the accompanying metal absorption systems 
amounting to a typical galactic rotational speed. There have been a few reports on the 
detection of the Lya emissions in the spectra of DLAs (e.g. Djorgovski et al. 1996, 1997; 
Moller & Warren 1993; Warren & Moller 1996). 

The first firm candidate is DLA 2233+131 reported by Djorgovski et al. (1996). The 
emission feature is characterized by a sharp blue edge, an absorption trough in the blue 
part, and a small secondary peak at around 5040 A. Djorgovski et al. (1996) ascribe the 
absorption feature to an ambient gas. 

On the other hand, Lu et al. (1997) suggest that some Lya forest clouds at lower 
redshifts may be responsible for the absorption feature in the blue part. However, as Djor- 
govski et al. (1996) and Lu et al. (1997) point out, it is highly probable that DLA 2233+131 
is a normal and isolated spiral galaxy. Therefore it will be a very rare coincidence if we 
have a Lya emitting DLA system obscured by other intervening clouds at a smaller red- 
shift. It is more probable, therefore, that the absorption feature is formed by the material 
surrounding the emission source. Lu et al. (1997) also comment that the absorption by 
either dust or the resonance scattering by neutral hydrogens in the galaxy can significantly 
alter the energy distribution of the escaping Lya photons. 

It is known that Lya photons from a galaxy are mainly emitted by the recombination 
processes in H II regions. Out of 8 nearby H II galaxies, Kunth et al. (1996) find 4 galaxies 
exhibiting a P-Cygni profile in the Lya emission. Subsequently, for one of such blue 
compact dwarf galaxies, Haro 2, Legrand et al. (1997) compared the observed P-Cygni 
Lya profile with the profiles properly scaled from the observed Ha profile of Haro 2. In 
their study it was found that the Lya profile shows a significant excess in the red part. 
They interpreted the red excess as the frequency redistributed flux mainly caused by the 
line photons back-scattered from an expanding shell. They conclude that the galaxies with 
a P-Cygni Lya profile may possess superbubbles. 

Furthermore, P-Cygni profiles are detected even at 3 < z < 3.5 in the normal star- 
forming galaxies (Steidel et al. 1996b). The similar P-Cygni profiles are also detected in 
the spectra of galaxies found in the Hubble Deep Field, where the redshifts are estimated 
to be z ~ 3 (Lowenthal et al. 1997). 

Therefore, we may draw a tentative conclusion that the P-Cygni type absorption 
feature appearing in the Lya emission of DLA 2233+131 strongly implies that the DLA 
system also contains similar superbubbles, which are often found in nearby blue compact 
H II galaxies. 

From the width of the absorption trough we estimate the bulk velocity of the expanding 
shell in DLA 2233+131 to be AV bulk « 200 km s _1 (see Fig. 2 in Djorgovski et al. 1996). It 
is observationally challenging to measure the local thermal velocity vth in the expanding H I 
region. The local velocity should also include the microturbulent motion, which may not be 
constrained observationally. Therefore, we assume that vth is in the range of 1 — 20 km s _1 
and adopt v t h ~ 10 km s _1 in this study. 

The slope at the blue edge of the primary peak discussed in the previous section is 
computed in the case of DLA 2233+131 observed by Djorgovski et al. (1996). According 
to them, the observed Lya line flux F\2iQ is 

F1216 = (6.4 ± 1.2) x 10" 17 erg cm" 2 s _1 . (26) 
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and the bin size is given by AA ~ 0.8 A. The maximum specific flux is F v ~ 8.5 juJy, 
which is converted to F\ ~ 1.0 x 10~ 17 erg cm -2 s _1 A . The number of bins between 
the blue edge to the maximum peak point is 6. Using our normalization scheme the slope 
is computed as 

10 4 F A AA 

'2233+131 = ^ -Z-T-. ~ ^50 (27) 

F1216 6 bins 

This result highly implies that the Sobolev optical depth To for the outflowing H I in 
DLA 2233+131 is a few times 10 5 , which is converted to the column density N HI ~ 
10 20 cm -2 assuming that s max does not greatly exceed 50. From these values and the 
overall shape of the profile, we exclude the possibility that Njji > 5 x 10 20 cm -2 or 
r >10 6 . 

Assuming a standard Friedman cosmology with Hq = 75 km s _1 Mpc -1 and fio = 0.2, 
the angular diameter distance to the DLA is 1.96 x 10 28 cm, and the luminosity distance 
tl = 8.14 x 10 28 cm. We suppose that the H II region is a sphere of radius R. Then the 
observed Lyct line flux is calculated by 

Fau = NMgj Vmu (28) 

r L 

where n e is the electron density, n p is the proton density, Vhii is the volume of the H II 
region, and the case B recombination coefficient cxb = 2.59 x 10 -13 cm 3 s _1 for T = 10 4 K 
(Osterbrock 1989). From Eqs. 26, 28, we obtain the electron density n e , 

1 k P c V _.__-3 



n e ~l [—j^J cm ( 29 ) 

We extrapolate the correlation between velocity dispersion and radius of giant H II regions 
(Terlevich et al. 1981) and use the velocity dispersion ~ 300 km s _1 of the H II region in 
DLA 2233+131 to estimate its size R ~ 1 kpc. If we assume that the emission region is 
fully ionized, then we get n e ~ 1 cm -3 , which leads to the total mass of the shell 

M HII ~ 1O 8 M . (30) 

Because one 05 star typically generates UV photons at a rate of 5 x 10 49 photons s -1 
(Spitzer 1978), it follows that the number of 05 stars needed to account for the Lya flux 
in DLA 2233+131 is ~ 10 3 under the assumption that the H II region is ionized purely by 
these 05 stars. According to Kennicutt et al. (1989) the first-ranked H II regions found in 
nearby spirals and irregulars require the total mass of the ionizing stars (> 1OM ) ranging 
10 2 -10 6 M Q . From this, we may conclude that the Lyct emission line from DLA 2233+131 
is originated from these first-ranked H II regions. 

4. Summary and Discussion 

We studied the line transfer in an optically thick and moving media including the 
damping wing scattering using a Monte Carlo method. The emergent profiles are charac- 
terized by the asymmetric double peaks with a primary peak in the red part, a secondary 
peak in the blue part, and an absorption trough in the blue. 
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The underlying mechanisms for the transfer and the dependence of the emergent pro- 
files on the H I column density and on the kinematics were investigated. The primary red 
peak recedes to the red as Nhi increases and for thicker media with To > 10 6 , the Lya line 
profile becomes two separated broad humps. It is also found that as s max gets smaller, 
symmetric double peak profiles are obtained, which is regarded as the static limit where 
bulk motion is negligible. We briefly discuss how to constrain the physical parameters tq 
and s m ax from an observed P-Cygni type profile. 

We have applied this method to explain the P-Cygni type profile in the Lya emission 
of DLA 2233+131 assuming that it is caused by an expanding medium with a high H I 
column density. The slope at the blue edge of the primary peak gives a strong constraint 
on r and we find that r ~ 10 5 . Noting that the expanding bulk velocity AVb u ik ~ 
200 km s _1 , the column density responsible for the P-Cygni type feature is deduced to 
be Nhi — 10 20 cm -2 . We also present physical quantities regarding the H II region and 
deduced that n e ~ 1 cm -3 , R ~ 1 kpc, and Mun — 10 8 M©, assuming near complete 
ionization. About 10 3 massive young stars are need to account for the Lya line flux of 
DLA 2233+131, is comparable with the first-ranked H II regions found in nearby spirals 
and irregulars. 

We briefly review the physical quantities of several astronomical objects which may 
possess similar kinematic properties to those of DLA 2233+131 discussed above. 

The first example we consider is the H I supershells and the worms/chimneys in our 
Galaxy (Heiles 1979, 1984: Koo et al 1991, 1992) and nearby galaxies. An especially large 
H I supershell is found in the nearby normal spiral galaxy M101 (Kamphuis et al. 1991). In 
M101 the measured line-of-sight expanding velocity component of the supershell amounts 
to ~ 50 km s _1 , which is smaller than that of DLA 2233+131 by a factor of ~ 5. The H I 
column density of the supershell is estimated to be N HI ~ 2 x 10 20 cm -2 . 

The second candidate is blue compact dwarf galaxies showing similar P-Cygni type 
profiles in their Lya emission (Legrand et al. 1997, Lequeux et al. 1995, Kunth et al. 1996, 
Puche et al. 1992). Kunth et al. (1996) found that 4 out of 8 observed H II galaxies exhibit 
P-Cygni profiles. One of them is Haro 2 whose H I column density is 2 x 10 20 cm -2 and the 
expanding bulk velocity ~ 100 km s _1 (Legrand et al. 1997). As Legrand et al. suggested, 
the possibility for the H II galaxy to exhibit the Lya emission or a broad absorption trough 
in their spectra depends on geometrical properties of the expanding shell, which include 
the angle between line-of-sight and the symmetry axis of the shell. It is noticeable that 
these blue compact galaxies are under active star formation processes. 

Currently, there are a few observational programs to discover primeval galaxies at 
z > 3. 

One is to search for star-forming galaxies like those found in 3 < z < 3.5 by Ly- 
man break characteristics (Steidel et al. 1996a, 1996b, 1997; Pettini et al 1997). The 
young galaxies detected in this way exhibit the similar P-Cygni type profiles to that 
found in DLA 2233+131. According to a deeper and higher resolution imaging (Giavalisco 
et al. 1996), these objects appear to be compact spheroidal cores often surrounded by 
low surface brightness nebulosities. Assuming the standard cosmology of O = 1 and 
H = 50 km s -1 Mpc -1 , a typical size of the cores is a few kpc. The authors suggested 
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that these are bulges of the primeval galaxies which are precursors of nearby normal galax- 
ies. 

Another program is provided by gravitational lenses. This method is now becomes 
a powerful tool, because the magnification of gravitational lens enables us to detect more 
distant objects than other method does. For example, the similar objects showing P-Cygni 
type profiles are also found in the Hubble Deep Field at z ~ 3 (Lowenthal et al. 1997). 
Trager et al. (1997) and Frye et al. (1997) also found a very distant objects at z ~ 4. 

It is remarkable that these objects show P-Cygni type Lyct profiles. Especially the 
galaxy discovered by Frye et al. (1997) shows several bright spots in the lensed image. A 
particularly notable point is that only one of them (i.e. N4 spot) shows a P-Cygni type 
Lya line profile whereas others show damped Lya absorption profiles (Bunker et al. 1997). 
This may tell us that there are several well-localized H II regions in the galaxy and only 
part of them shows a P-Cygni type Lyct emission, while the Lya photons of others are 
shielded by an optically thick medium along the line of sight. Analogous situations are 
provided by nearby starburst dwarf galaxies or supershells, and we suggest that these 
starburst dwarf galaxies are the very counterparts of the above three kinds of primeval 
galaxy candidates including the DLA candidates. 

Such outflowing media seen in extragalactic objects are often accompanied by the 
sites associated with active star formation (Reach et al. 1993). Because DLA 2233+131 
and other remote galaxies are believed to be prototype of present galaxies, it is very 
interesting that there may exist primeval galaxies at 3 < z < 5 showing starbursts, which 
also gives constraints on the cosmological models for the structure formation. Moreover 
this speculation supports the picture that star formations are continuous processes from 
z = to z ~ 5. 

We propose complementary observations of metal lines in the spectrum of DLA 
2233+131 in order to provide independent constraints on the relevant physical quantities. 
We also propose the further imaging of DLA 2233+131 by the Hubble Space Telescope 
which may provide independent and important clues to the identity of the system. 

In a subsequent paper, we will investigate in more detail the detectability of the 
Lya emission in these star forming objects, which is thought to be related with the geo- 
metrical factors such as inclination associated with the shielding of the dust lane. 

In numerical approach it will be desirable to develop a more flexible code to treat 
inhomogeneous media with an arbitrary velocity field and various shapes to secure the 
constraints on the characteristic physical quantities such as the column density and the 
bulk flow scale. 
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FIGURE CAPTION 



Fig. 1 

A surface plot of the s(x,t) table, x = Av/Avd and TAU = r represents the optical 
depth. The plateau on the right top side corresponds to s = 2s max + 1. 

Fig. 2 

The emergent Lya profile from a thick expanding medium. The horizontal axis represents 
A A/ A Ad and the vertical axis for the relative flux. The emission source is assumed to be 
given by a Gaussian profile oc e~^ x ^ 2<Tx ' , where the width a x is set to be 5. We choose 
Smax = 10 and tq = 10 4 . The emergent profile is represented by the solid line and the 
dotted line shows the initial Gaussian profile. 

Fig. 3 

The escaping optical depth and the single-scattering escape probability. The dashed line 
represents the optical depth r x , and the escape probability P esc = e~ Tx is shown by the 
dotted line. The dot-dashed line is the input emission profile at the origin given by a 
Gaussian with the width a x = 10. The solid line is the product of the initial emission at 
the origin and the escape probability. 

Fig. 4 

The emergent profiles for a flat incident continuum. The parameters are the same as in 
Fig. 2. 

Fig. 5a 

A typical scattering behavior of a photon with x ~ —5 The horizontal axis stands for 
the scattering number. In the bottom frequency deviation is shown, and in the top panel 
the scattering type is given, where stands for a wing scattering and 1 for a resonant 
scattering. In the middle panel the radial distance of the photon from the center of the 
scattering sphere is shown. 

Fig. 5b 

Another typical scattering behavior of a photon with x ~ 5. See the text for more detail. 
Fig. 6 

The emergent Lya profiles for various bulk velocity scales s max of the scattering medium. 
We fix ro = 10 4 and the initial Gaussian profile having a width of a x = 5 is used. The 
thick solid line is for the case of s max = 0.5, the dotted line for s max = 2, the long dashed 
line for s max = 1, the dashed line for s max = 5, and the dot-dashed line for s max = 10. 

Fig. 7a 

The emergent Lya profiles for various column densities of the scattering medium. s max = 
10 and the initial Gaussian profile having a width of a x = 5 is used. The thick solid 
line indicates the emergent Lya profile for r = 10 4 which corresponds to N HI = 4.1 x 
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10 18 cm 2 , the dotted line for r = 10 3 , and the dashed line for r = 10 2 . The small dips in 
the secondary peaks are numerical artifacts caused by the truncated matter distribution. 

Fig. 7b 

The emergent Lya profiles for various column densities of the scattering medium. The 
thick solid line indicates the emergent Lya profile for To = 10 4 which corresponds to 
N m = 4.1 x 10 18 cm -2 , the dotted line for r = 10 5 , and the dashed line for r = 10 6 . 
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